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We show that the energetic instability of dark solitons 
is associated with particle-like motion, and present a simple 
equation of motion, based on the Mossbauer effect, for dark 
solitons propagating in inhomogeneous Thomas- Fermi clouds. 
Numerical simulations support our theory. We discuss some 
experimental approaches. 

PACS number(s): 03.75.Fi, 03.65 Ge 



The success of the Gross-Pitaevski mean field theory 
in describing experimentally observed dilute Bose con- 
densates shows that one really can persuade a large 
number of particles to behave as a field. The perversely 
obvious thing to do with such a field, then, is to persuade 
it to behave as a particle. This would not only be a pleas- 
ant closing of a circle, but would link condensate physics 
to other subjects in which solitons arc important, such as 
fluid mechanics, nonlinear optics, and fundamental par- 
ticle theory. Perhaps most importantly, it might offer a 
powerful laboratory for studying the interface between 
quantum and classical physics, by providing a wholly ar- 
tificial classical particle, composed of highly controllable, 
weakly interacting quantum particles. In this Letter we 
discuss one particular particle-like configuration of the 
Gross-Pitaevski mean field, namely the one-dimensional 
dark soliton. We describe the behaviour of dark solitons 
in inhomogeneous potentials, reporting both analytical 
and numerical results, and we propose some experimen- 
tal approaches for effectively one-dimensional traps. We 
also suggest that the dark soliton should be considered 
the simplest member of a family of technically unstable 
but nevertheless robust monopoles in D dimensions. 

The Gross-Pitaevski equation (GPE) governs the evo- 
lution of the c-number 'macroscopic wave function' 
^{x^t) of a Bose-Einstein condensate. By appropriately 
scaling the wave function, space, and time, and incorpo- 
rating a chemical potential by extracting a factor e"*** *, 
one can write this equation in the convenient form 



idtijj 



-vV^ + (IV-P + y(x) - k2)V' 



(1) 



We have assumed a positive scattering length, and have 
not restricted the normalization constant Jrixji/'P- We 
will focus for most of this paper on the limit of a long trap 
so thin that one can apply the GPE in one dimension. 
The approach to this limit from three dimensions has 
recently been discussed 

Eqn. ([^) in one dimension with V = has been ex- 
tensively studied in nonlinear optics and a particle- 
like solution has long been known: ^ds — Ktanh(/ta:). 



This time-independent solution is known as a dark soli- 
ton, because it describes a small dark spot in a light 
pulse. If '4!{x) were restricted to be real, the dark soli- 
ton would be topologically stable in the same way as the 
Sine-Gordon kink that it strongly resembles. Since ip is 
very nearly constant far from x = 0, but changes sign 
within a few healing lengths, it resembles a domain wall: 
a one-dimensional monopole. 

There exist close relatives of the dark soliton which 
may be recognized as D-dimensional monopoles; we di- 
gress briefly to sketch the family group. If we allow ip 
to be a multi-component object ip = {ipi, 4'd)i and in- 
terpret IV'P = J2j IV'jP: then (0) is the mean field equa- 
tion of motion for a condensate of atoms with D internal 
states (among which all scattering lengths are equal). In 
two dimensions one finds "0 = /(nr) {cos 6, sin 6) to be a 
solution to the two component GPE, where / is the ra- 
dial profile of a one-component vortex. Analytical and 
numerical results related to this solution will be reported 
elsewhere. The generalization to three dimensions for 
three components is obvious In all these cases, if ip 
were restricted to be real we would have topologically 
stable monopoles, as small 'bubbles' in the condensate. 

It is not difficult to see, however, that by taking tp 
into the complex plane one can continuously deform it 
into a configuration with constant density, eliminating 
the 'bubble' in One can exhibit this instability 

more precisely by exploring the free energy of nearby 
configurations; although this analysis can be carried out 
more generally, we now revert to one dimension. As a 
time-independent solution to (|^) {V = 0, one dimension), 
ipDS = tanh Kx is a stationary point of the free energy 
functional 



G 



dx [l^f + (|Vp-«')' 



(2) 



One can diagonalize the Hessian matrix of G at the sta- 
tionary point by solving what may be recognized as imag- 
inary time Bogoliubov equations for the perturbation 6ip. 
(The real time Bogoliubov theory also diagonalizes the 
Hessian matrix, but by a symplectic instead of orthogo- 
nal transformation, as required by canonicity.) Because 
-0 is real, the real and imaginary parts of 6ijj = R + iS de- 
couple, and independently satisfy (different) Schrodinger 
equations: the normal modes Rk{x) and Sk{x) are re- 
spectively eigenstates of and Hi, for 
1 d2 



2dx^ 



nip 



2 

DS 



(3) 



One finds the 'core bound state' 5*0 = 1/ cosh^a;, for 
which the Hessian matrix has a negative eigenvalue; this 
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shows that there are configurations near the dark soliton 
with lower free energy. And in numerical evolution of the 
dark soliton in imaginary time, numerical noise soon finds 
this unstable mode, and fills in the monopole bubble. 

Imaginary time evolution has only the most indirect 
implications for real time evolution, however. In one di- 
mension with V = constant, we have the following 
two parameter family of exact, but not generally time- 
independent, solutions to (|l|): 

ipq,q{x, t) = iq+ ^ - (p' tanh \J n"^ — q'^{x - q) , (4) 

where q is the time derivative of q, constant and satisfying 
satisfying \q\ < k. 

For these moving monopoles (familiar in nonlinear op- 
tics as 'grey solitons'), IV'P never drops below q^; and the 
phase slip across the soliton is tt — arctand^jZ-^/K^ — ip). 
This means that in the limit q ±k, where the soliton is 
travelling at the superfluid critical velocity, the dark soli- 
ton actually vanishes, becoming identical with motionless 
condensate. But the free energy of the grey soliton is 

G{q,q) = ^{K^-e)^ . (5) 

This implies that grey solitons have negative kinetic en- 
ergy. The imaginary time Bogoliubov equations found 
only one negative energy degree of freedom near the dark 
soliton, and we can now identify it in the real time Bo- 
goliubov theory as the canonical momentum conjugate 
to the translational zero mode. And so we see that the 
instability of the dark soliton is not merely to break up or 
fill in, but to acquire velocity. 

Since dark solitons have both translational modes and 
independent velocities, they have the phase space of a 
particle. This is in contrast to a vortex in two dimen- 
sions, whose two positional degrees of freedom are in fact 
canonically conjugate, so that the vortex phase space is 
only two dimensional |^]. It is therefore natural to seek 
a particle-like second order equation of motion for dark 
solitons in an inhomogeneous potential V{x). One can 
readily guess that, at least in a sufficiently slowly vary- 
ing potential, one should simply replace — ^ — V{q) 
in (^, and insert the result into a background Thomas- 
Fermi cloud. The question of how such a configura- 
tion will evolve in time, however, is not trivial. Time- 
dependent variational approaches can be informative but 
also misleading. Fortunately, a more systematic proce- 
dure exists. 

In the experimentally relevant situation where V{x) 
varies slowly on the healing length scale one can 

apply a time-dependent boundary layer theory similar 
to that used for vortices in Ref. [Q. Perturb around 
the Thomas-Fermi approximation everywhere outside a 
small moving 'core zone' extending a few healing lengths 
from the monopole; treat the inhomogeneity of 1^ as a 
perturbation within the core zone; and then match solu- 
tions smoothly at the moving boundary. This procedure 
leads to a very simple result in the case where the only 



excitations present (other than the soliton motion) are 
those generated by the soliton as it accelerates or de- 
celerates. In this case the boundary layer theory demon- 
strates that when the healing length is much smaller than 
the Thomas-Fermi cloud, the system is in the Mossbauer 
limit, in which phonon mode excitations conserve mo- 
mentum without contributing significantly to the total 
energy. As a result, the soliton energy [k^ — q^ — V{q)]'^/^ 
is conserved to leading order in the ratio of length scales, 
and this implies the equation of motion ||] 

Q^-lv'{q). (6) 

In a harmonic trap, Eqn. (^ implies oscillation of the 
soliton 1^ with frequency 1 / a/2 times that of the dipole 
mode of the condensate. We have confirmed the fre- 
quency factor to rather more than the expected accuracy 
in numerical simulations of harmonic traps over a wide 
range of condensate densities and oscillation amplitudes; 
we have also monitored the center of mass motion, and 
confirmed that it is decoupled and has the trap frequency. 
Eqn. (^) also holds for arbitrary potentials, however, as 
long as they vary slowly on the healing length scale. We 
have therefore further confirmed the good accuracy of our 
Mossbauer limit equation of motion by solving Eqn. (|l|) 
numerically over a wide range of parameters and for var- 
ious potentials; a generic sample is shown in Fig. 1. 

To solve the NLSE numerically we use the split-step 
Fourier method [Q, in two stages. First we construct an 
initial state with no excitations apart from the soliton, 
by propagating a crude trial wavefunction in imaginary 
time, and re-normalizing it between each step. Eventu- 
ally, numerical noise excites the 'core bound state' and 
the monopole rapidly disappears in imaginary time; but 
before this there is a clear plateau in energy, during which 
the monopole slowly drifts in the Thomas-Fermi cloud, 
in the direction of decreasing density. We simply stop the 
imaginary time evolution during this period, and save the 
configuration. Then we propagate this optimized initial 
state in real time. As is well known, real time evolution is 
prone to (clearly recognizable) high frequency numerical 
instabilities; but by keeping our time steps shorter than 
/i^/tt, where h is the spatial grid interval, we are able 
to simulate several periods of monopole motion, even at 
condensate densities approaching those attained in real 
experiments. (At /dxjV'P = 1000, our Thomas-Fermi 
cloud extends to over ten times the trap width.) 

Since Eqn. (^ predicts a particular frequency for fi- 
nite amplitude oscillations around a local minimum of 
the trap potential, and a dynamical instability for mo- 
tion near a local maximum, it would be interesting to 
compare these predictions with those from Bogoliubov 
theory. Since we have numerically found frequencies very 
close to that from Eqn. (|^) even for extremely small am- 
plitude oscillations, we strongly suspect that the two the- 
ories agree in the Mossbauer limit. Such agreement is not 
actually necessary, however: Eqn. (|6|) is derived for ar- 
bitrary amplitude oscillations but neglecting coupling to 
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phonon modes, while the Bogohubov spectrum governs 
only infinitesimal perturbations, but includes all excita- 
tions. Co-incidence of the two spectra would imply that 
the dark soliton decouples from phonon modes at ex- 
trema of the potential, a result which might be of some 
importance for the problems of dark soliton zero point 
motion and mesoscopic tunneling. 
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FIG. 1. Normalized {]dx\il}f = 300) for a 

dark soliton oscillating in an asymmetric potential with 
a bump, shown (magnified) in dots. The potential is 
V — 0.1a;(a; — 2) -|- 1.1/ cosh^(a:), the soliton starting point is 
x=-2.81 (with a grid uncertainty of about 0.03); Eqn. (|^ pre- 
dicts the interval between successive times of complete dark- 
ness (IV'P = at the minimum) to be r/2 = 11.5 ± 0.1. 
Comparing times in lower left and lower right plots shows the 
Mossbauer limit is good to within 5%; the ratio of healing 
length l/l^l to potential scale is of this order. 

Since the equation of motion (|^) is obtained by ne- 
glecting excitation of phonon modes by the monopole 
motion, but in fact some excitation does occur, we must 
expect dissipation (which could in principle be calculated 
by continuing the time-dependent boundary layer theory 
to higher orders in the ratio of length scales). Since the 
soliton energy is negative, however, dissipation has the 
effect of increasing the amplitude of the monopole mo- 
tion in the trap. We have observed this effect in numer- 
ical simulations: phonon mode excitations are visible in 
slight distortions of the Thomas-Fermi envelope. These 
excitations produce small variations, plus a slow trend 
upwards, in -|- ^(g); see Fig. 2. 
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FIG. 2. Plots of the quantity -\- F(q), conserved in the 
Mossbauer limit, versus time, as a dark soliton moves in the 
potential of Fig. 1. The three curves are for various densities: 
\dx\'4if = 100,500,1000 from top to bottom. For we 
actually take the height of at the soliton minimum, which 
is equal to for dark solitons in bulk. 

We conclude our theoretical discusion by remarking 
on the interaction of two solitons. Analytical work on 
dark solitons in bulk has shown that two dark solitons 
interact through a short range repulsive potential, whose 
maximum height is finite and velocity dependent [ pd] |. 
So solitons more than a few healing lengths apart do 
not influence one another, and solitons colliding rapidly 
can pass through each other with negligible interaction; 
but two slow solitons cannot pass each other. Moreover, 
the transition between impenetrability and near-freedom 
is fairly sharp as a function of relative velocity. On the 
other hand, the presence of a second soliton in a harmonic 
trap appears to shorten the oscillation period of the pair; 
see Fig. 3. 
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FIG. 3. Two difi'erent cases of two solitons oppositely dis- 
placed in a harmonic trap. The trap angular frequency is 
one, and the time axis is shown in units of 2n. The solid line 
shows solitons passing through each other; the dashed line, for 
two solitons initially displaced somewhat less, shows bounc- 
ing motion. For both cases, Jdx\ip\'^ = 100. Note that the 
period of bouncing motion is approximately 2-k x 1.33, that 
of passing motion 2n x 1.36. For comparison, the period of a 
single soliton at this density is 2-k x 1.42. Both motions are 
also clearly modulated, due to coupling to a phonon mode. 

We now briefly outline experimental approaches to the 
creation of dark solitons. An ingenious proposal has re- 
cently been presented to construct dark solitons by adi- 
abatic passage This scheme requires that the con- 

structed soliton be a dynamically stable stationary state, 
which would seem to prevent one from making dark soli- 
tons on slopes or local maxima of the potential. Since the 
capability exists to make temporary small local wells any- 
where within the larger trap, however, this need not ac- 
tually pose a problem. (This capability should also allow 
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more general testing of the dark soliton equation of mo- 
tion, as one could arrange a sequence of 'speed bumps', 
or accelerating wells, along the thin trap.) Other recent 
work has noted that dark solitons are generically pro- 
duced during collisions between condensates |l^Jl^, or 
between generic travelling disturbances . 

Another challenge is to detect dark solitons once they 
are made, since they are typically much smaller than a 
wavelength of light. Current techniques of relaxing the 
trap so that the soliton 'hole' grows to visibility should 
be helpful, as should the interferometric result that a 
dark soliton produces a discontinuity in the pattern of 
Ramsay fringes. We would like to suggest an alternative 
method, however, which could allow one to view dark 
solitons directly and non-destructively: trapping atoms 
with trapped atoms. Since there is a bound state in the 
dark soliton potential, a few atoms of another species 
could in principle be bound inside the dark soliton as 
it moves (assuming that they are repelled by the con- 
densate atoms). In principle one might hope to intro- 
duce these probe atoms simply by 'dusting' the conden- 
sate with many of them, so that most probe atoms are 
simply repelled into the surrounding thermal cloud, but 
some are sympathetically cooled and 'stick' inside any 
dark solitons which may exist. Once there, they may 
be made to fluoresce independently of the condensate, so 
that the 'dark' solitons may be followed as small bright 
spots. Further investigation of this proposal is obvi- 
ously required, as some condensate atoms are lost as the 
probe atoms are cooled sympathetically, and one must 
ensure that the capture cross section of the soliton is 
high enough that one does not need so much 'dust' that 
the soliton, or even the condensate, is destroyed. In prin- 
ciple however it seems a promising way to observe many 
types of topological defects. 

As a final experimental remark, we note that in traps 
which are not sufficiently close to the one-dimensional 
limit one expects dark solitons to be unstable, presum- 
ably above some threshold velocity, to the formation of 
vortex rings; these would of course be at least as inter- 
esting as pure solitons. 

In conclusion we wish to re-emphasize that although 
dark solitons are unstable, this instability should be much 
more of an advantage than a problem, as it allows the 
solitons to behave as particles, with independent posi- 
tion and velocity. In the Mossbauer limit a very simple 
equation of motion applies, which is expected to remain 
valid while the amplitude of soliton motion in the poten- 
tial slowly grows. The final victory of the instability is 
the escape of the dark soliton from the trap: at least for 
these mesoscopic 'hollow men', the world will after all 
end with more of a bang than a whimper . 
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